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Abstract 

In this paper, we study the supercurrent in a mesoscopic Josephson junction 
(MJJ) and its quantum statistical properties in the presence of nonclassical light 
fields. We investigate in detail the influence of external nonclassical light fields 
on current- voltage step structures of the MJJ. We also study in detail quantum 
statistical properties of the supercurrent when the external quantum electro- 
magnetic fields are even and odd coherent-state light fields. It is shown that the 
supercurrent in the MJJ exhibits both squeezing effect and quantum coherences. 
It is demonstrated that the MJJ can feel the difference not only between classi- 
cal light fields and nonclassical light fields but also between different nonclassical 
light fields. 



1 Introduction 



In last few years considerable attention has been devoted to mesoscopic Josephson 
junctions (MJJ) [1-6] due to the prediction of new quantum effects [1,7-10] in these 
small-capacitance junctions. Much work concentrated on the macroscopic quantum 
tunnelling of the phase of current-biased junctions and the influence of dissipation on 
them [11-15]. It is well known that the effect of classical electromagnetic fields(EMF) 
on macroscopic Josephson junctions has been well understood. The Shapiro steps [16] 
are the result that the external static EMF and time-varying EMF work together on 
Josephson junctions. It is natural to consider the effect of quantum EMF on Josephson 
junctions. In particular, it is an interesting topic to study the effect of quantum EMF 
on the MJJ since it can respond to the quantum nature of the external EMF. 

On the other hand, nonclassical light fields [171 have been studied extensively in 
the recent years and produced experimentally. These studies emphasize nonclassical 
effects of light fields such as squeezing effect and sub-Poissonian distribution to reflect 
the quantum nature of light which can not be understood classically. Even and odd 
coherent states(CS)|ll8|, p!9| are two types of important nonclassical states of light and 
have wide applications in physics [21-24]. Although even and odd CS can be decom- 
posed into superpositions of the Glauber CS, their properties are quite different from 
those of the Glauber CS since the former are nonclassical states while the latter are 



classical ones|2^. On the other hand, even CS and odd CS exhibit significantly different 
nonclassical characteristic separately |1^, |2^ . Even CS exhibits quadrature squeezing 
but does not sub-Poissonian photon statistics while odd CS exhibits sub-Poissonian 
but does not quadrature squeezing. 

More recently, Vourdas [^ has studied the effect of external quantum EMF on the 
MJJ by the use of a circular superconducting device with a MJJ in the presence of 
quantum EMF. He has investigated the step structures of the current and voltage of 
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the MJJ. It has been shown that these step structures are similar to Shapiro steps in 
macroscopic Josephson junctions, but their details depend on the nature of quantum 
EMF. 

The aim of this paper is to study systematically properties of the supercurrent 
through the MJJ in the presence of nonclassical EMF by taking even and odd CS as an 
example of nonclassical EMF and discuss effects of the external nonclassical EMF on 
these properties concretely. We shall study the current- volt age step structures of the 
MJJ and quantum statistical properties of the supercurrent. This paper is organized 
as follows: in Sec. 2 the current-voltage step structures (CVSS) of supercurrents in the 
MJJ and the effect of nonclassical EMF on them are investigated. Sec. 3 and Sec. 4 
are devoted to study quantum statistical properties of the supercurrent. Squeezing 
properties and quantum coherences are investigated in detail for even and odd CS 
EMF. The results are summarized and conclusions are drawn in the last section. 

2 CVSS of the MJJ in the presence of even and 
odd coherent-state EMF 



Consider a MJJ with a capacitance C in a circular superconducting device |^. The 
MJJ is described by the following Hamiltonian {h = c = 1) 

H=^^±^ + Ej{l-cose) (1) 

where the first term is the charging energy, and the second one is the Josephson coupling 
energy; q is the charge on one of the electrodes of the MJJ, and Q is the external charge; 
6 is the phase difference of the superconducting order parameters in the two electrodes; 
and Ej is the Josephson coupling constant. 

Capacitance C of conventional Josephson junctions and the temperatures where we 
usually operate allow us to neglect the charging energy in (|I|) and to forget the quantum 
mechanical origin of the Josephson effect. In stead, we view the Josephson junction as 



a classical object whose state is characterized by a sharp value of the phase difference 
across the junction. However, In the case of the MJJ, the classical descriptions is no 
longer sufficient, and the charging energy is not negligible. We have to describe the 
MJJ in a fully quantum mechanical way. In this case, the charge q and the phase 
difference 9 are a pair of quantum mechanical conjugation variables with q = —i2ede, 
they satisfy the canonical commutation relation: 

[e,q]=z2e (2) 

The Hamiltonian (|l|) is then rewritten as 

H = ~E,ide + i^f + Ej{l - cosO) (3) 

where the Josephson coupling energy related to the critical current Icr and the charging 
energy associated with the Cooper-pair charge 2e are given by, respectively. 



Ej-Ye^ ^c-^ (4) 



Starting with the Hamiltonian (^, one can find that the Josephson-current operator 
/ in the Heisenberg picture takes the form: 

/(t) = /crsin^(t) (5) 

As is known, the Shapiro steps on macroscopic Josephson junctions is the result 
which the classical static electromotive force work together with the classical time- 
varying electromotive force. In the present case, we still impose the classical static 
electromotive force Vq but the quantum time-varying electromotive force on the MJJ 
through applying a classical magnetic flux (f) = VqI and a ac quantum single-mode 
magnetic field of angular frequency uji, respectively. Then, the external voltage induced 
by the quantum magnetic field in the Heisenberg picture can be written as 

K.(t) = ^(e*^^*a+-e--*a) (6) 



where a and a"*" are the annihilation and creation operators of the quantum EMF. 



we consider the case ooi ^ y/E^Ej where y/EcEj is the angular frequency of the 
MJJ when it is approximated as a harmonic oscillator. In this case the Hamiltonian 
of the MJJ itself produces an evolution which is much slower than the Hamiltonian 
describing the effect from external fields, and the phase evolves approximately in the 
form [27]: 

where Uq = 2eVo, ^ = \f2e (e is the charge of the electron), and -D(Z) is the displace- 
ment operator defined by 

D{Z) = e^"+-^*« (8) 

Then the expectation value of the Josephson current / can be given by tracing of the 
current operator / over the density operator of the external quantum EMF, 

I = /c,Im {e^'"°*Tr[pD(Ce^'"^*)]} (9) 

It is interesting to note that the density operator p of the external quantum EMF 
controls the operation of the MJJ in the sense that the Josephson current and also all 
the other quantities, one might wish to calculate, depend on the density operator p. 
So that we can change the state in which the junction operates through changing the 
density matrix of the external quantum EMF. This gives a way to realize the controlling 
of the Josephson tunnelling. 

We now consider the current- voltage step structures of the MJJ when the external 
quantum EMF are in even and odd coherent states, respectively, and investigate the 
effect of quantum nature of the external quantum EMF on the step structures. 

The number-state representations of even and odd CS are defined by 

oo y2n 

\Z), = Ar,(|Z|)$:^=|2n), N,i\Z\)^cosh-mZ\') (10) 



„=o J{2n)\ 



oo n'2n+l 

\Z), = NMZ\)Y1 I |2n+l), iV,(|Z|) = sinh-^(|Zp) (11) 



n=o y'(2n + 1)! 

which can be decomposed into superpositions of Glauber CS \Z) and | —Z) in the form 

= Ai{\Z \){\Z)+ I -Z)\ A-\\Z I) = 2(1 + e-^l^^), (12) 
= Ak\Z\){\Z)- \-Z)\ A-_\\Z\) = 2(1 - e-^l^^ (13) 

It should be noted that although even and odd CS have the above simple decomposi- 
tions of Glauber CS, they have quite different quantum nature since even and odd CS 
light fields are nonclassical light fields while the Glauber CS light field is the classical 
one 
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To get the Josephson current, one needs the following matrix elements of the dis- 
placement operator in the CS representation 

^±z,±z(ee^"^*) = {±Z\D{^e''-^')\±Z)=exp{-^e±^2u), (14) 

^±z,Tz(ee*"^*) = {±Z I D(ee^-^*) I tZ) = expi-^e - 2 |Z T ±2t^) (15) 

where 

u = \Z\ sin(ci;it — (p), v = ^ \Z\ cos(ti;it — ip). (16) 

Here we have taken Z =\Z\ e**^. 

Making use of Eqs. (p^ ) and (0), from the expressions (0) and (^ one can obtain 
the expectation values of the phase- difference operator in the even and odd CS with 
the results 

e{Z\e''^'^ \Z), = 2A+(|Z|)e'"«V5«'(cos2u + e-2^'cosh2t;), (17) 
o{Z\e'^^'^ \Z)o = 2A_(|Z|)e'"«V5«'(cos2n-e-2^'cosh2i;) (18) 

and the Josephson current in the presence of even and odd coherent-state EMF, 

le = 2IcrA+{\Z\)e-^2^\cos2u + e-^^\osh2v)smuJot, (19) 
lo = 2I^rA-{\Z\)e'^^\cos2u- e-^^\osh2v)smujot (20) 



In order to see the CVSS of the above currents, we expand the rhs of Eqs. (|T9[) and 
Of) as follows: 



oo 

.lt2 



h = 2hrA4\Z\)e-^^ J2 [J2n{2^\Z\) 

n=— oo 

+e-'l^'/2n(2e \Z\)] sin[(2ncui + uJo)t - 2nip], (21) 



oo 

lt2 



n=— oo 

-e"^'^'/2„(2^ \Z I)] sin[(2r2cji + cuo)t - 2rup] (22) 

where J„(x) and /„(x) are Bessel and modified Bessel functions, respectively, and we 
have used the formulae: 



e 



ixsino _ J„(a;)e'"^ e^^°^^ = ^ /n(x)e*"^ (23) 



The dc components of the Josephson currents can be obtained through taking the time 
meanvalues of the expressions (^TJ) and (^) in the interval (0,t). It can be seen from 
Eqs. (^) and (^) that for even CS and also for odd CS when the dc voltage Vq obeys 
the following condition 

2nuJi — ujq = 0, i.e., Vq = — (24) 

e 

we get Shapiro dc currents similar to those in macroscopic junctions. Their expressions 
take the form 

If = 2J„,A+(|Z|)e-^?^ [j_2„(2e \Z\) + e-2|^^/_2„(2e \Z\)] sin(2n<^), (25) 
/f = 2hrA4\Z\)e--^^' [j„2n(2e \Z\) - e-2|^'/_2n(2e 1^1)1 sin(2ny;). (26) 



It can be seen from Eqs. (p^f), (p5|) and (^) that even CS and odd CS have similar 
CVSS. Especially, they share the same voltage step structures when Lp = ^ (n ^ 
0, k = 0, ±1, ±2, . . .), while the dc components vanish, i.e., I^^ = /^^ = 0; and when 
|Z 1^ 1, we find that even CS and odd CS have almost equal dc components, that 
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is Jf = = Jc^e-^^V_2n(2^ \Z I) sm(2n^). Thus, when \Z |> 1, even and odd CS 
have the same CVSS. Nevertheless, when |^ |^ 1, from Eqs. (^) and (^) we have 
/f = and /f = i/cre-^«' [J-2n(2^ |Z|) + /_2„(2^ |Z |)] sin(2n(^) which indicate that 
when |Z 1^ 1 even and odd CS have quite different CVSS. Therefore, the MJJ can feel 
the difference between different external nonclassical EMF. 

On the other hand, the CVSS of even and odd CS are also different from that of the 



Glauber CS |27|. The voltage jumps (0) in the former are double in size in comparison 
with the voltage jumps in the latter. This means that the electrostatic energy that the 
Cooper pair loses during tunnelling in the case of even and odd CS is twice of that in 
the Glauber CS case. As concerns the dc components of the Josephson currents, the 
dc components for even and odd CS have only one direction for all steps for a given Z 
while for the Glauber CS the dc components have two directions which correspond to 
n being even and odd, respectively. Taking into account the results in Ref. [^], we can 
find that the Glauber CS, CS with thermal noise, and CS with partially randomized 
phase have the same CVSS, while even and odd CS, squeezed vacuum states exhibit 
the same CVSS. This reflects the fact that the former three types of states are classical 
states while the latter three ones are nonclassical states. Therefore, we can conclude 
that the MJJ can also feel the difference between the external nonclassical EMF and 
classical ones. 

It is worthwhile mentioning that the existence of the CVSS of the MJJ depends 
strongly on the form of the density matrix describing the external quantum EMF. 
Those states which off-diagonal elements of their density matrices in the number-state 
representation are nonzero, such as the Glauber CS, even and odd CS, squeezed vacuum 
states, exhibit the CVSS of the MJJ, while for those states which off-diagonal elements 
of the density matrices are zero, such as CS with randomized phase, thermal states and 
number states, there do not exist the CVSS. This indicates that the existence of the 
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CVSS of the MJJ are determined by the external dc voltage Vq and the off-diagonal 
elements of the density matrices of the external quantum EMF. In other words, the 
step structures depends on the dc voltage Vq and the quantum interference among the 
number states of the external quantum EMF. In this sense, the step structures reflect 
to some extend quantum coherences of the external quantum EMF. 

3 Squeezing properties of supercurrents in the MJJ 
in the presence of quantum EMF 



Squeezing is one of the most important concepts in quantum optics p5| , p8| . It has 
potential applications in low-noise optical communications and high-precision interfer- 
ometic measurements. It is of importance to study squeezing behaviors of phase in 
quantum mechanics p9| , |30| . As is known, the supercurrent in a Josephson junction is 
a quantum-phase- correlation phenomenon. So it is interesting to study squeezing prop- 
erties of the supercurrent. In this section, we discuss squeezing of the supercurrent in 
the MJJ in the presence of quantum EMF. 

Consider two noncommuting Hermitian operators A and B acting on the state 
space of the external quantum EMF. The variances (fluctuations) of these operators 
((Ai)2) = {A^) - and ((AB)^) = (S^) - {B)"^ obey the Heisenberg uncertainty 
relation 

((Aif)((AS)^)> (27) 

which is fundamentally different from that of the position and momentum operators 
of a quantum mechanical particle ((Ax)^) ((Ap)^) > ^ since the commutator [A, -B] 
on the rhs of equation(|27|) is not a C number but an operator which makes the rhs 
of the uncertainty relation (p^) to be state-dependent. Uncertainty relation with a 
state-dependent rhs is well known from earlier studies of the atomic coherent states 



Following Wodkiewicz and Eberly we shall say that the variances of 
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the operators A and B are squeezed if 



((Aif)<i|([i,fi])| or ((A5)^)<l|([i,E]) 



(28) 



In order to measure the degree of squeezing of these operators, we introduce two squeez- 
ing parameters 5*^ and 5*^ defined by 



((Ai)^)-||([Ag])| ((A5)2)-i|([i,fi]) 

— 1 , „ I , , 



1\{[AB])\ ' - II ([A 5]) 

The squeezing condition now takes the simple form: 
Sa<0 or < 



(29) 



(30) 



-1 (X 



The maximum(100%) squeezing of variance ((AX)^) corresponds to Sx - 
A or B). 

For the present case of the MJJ in the presence of quantum EMF, from the canonical 
commutation relation of the charge on the electrodes and the phase difference operator 
(^, we can obtain the commutation relation between the supercurrent operator and 
the charge operator, 



[/, q] = —i2eIcT cos 9 

Then, the degree of the squeezing of the Josephson current is given by 

((Asin^)2) -e |(cos^)| 

S[ — 

e I {cos 6) I 

For simplicity we have taken Jcr as the unit of measuring the supercurrent. 
For even and odd CS, it is easy to get 



(31) 



(32) 



3'"''* [^z,z(ee'"^*) + I^_z,z(ee'"^*: 
+ Dz,-z(ee^"^*) + /^-z,-z(ee^"^*) 
5*""* [/^z,z(ee*"^*) - /^-z,z(ee*"^*: 



(33) 



(34) 
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where the matrix elements of the displacement operator are given by Eqs. (|1^) and 
From Eqs. (BSf) and (Bl) we find that 



e{Z\ I ^IJ^^ j \Z)e = 2A4\Z\)e--^^' [cos2M + e-2|^'cosh2t; 



o{Z\ I '''''i ] \Z)o = 2A(\Z\)e-^^' fcos2M-e-2|^'cosh27; 
\ cos 6* / L 



sin Uot 
cos Uot 

sin uot 
cos Uot 



(35) 



(36) 



Then, we obtain the fiuctuations of the supercurrent in the MJJ for even and odd CS 



((Asin^)2)„ 



- — S+ cos 2uot — sin^ u^t, 

- — cos 2uL!Qt — sin^ uj^t 



(37) 
(38) 



where 



B± = A±{\Z\)e-^^\cosAu±e-^\^\osh4v), 
C± = 2A±{\Z\)e-^^^\cos2u±e"^^\osh2v). 

On substituting Eqs. (|35|)-(|38D into Eq. (^), we find the squeezing degree of the 
supercurrents with respect to even and odd CS, respectively 
1 - 2B+ cos 2uJot - 2Cl sin^ uot 



2e I C+ cos uot I 
1 - 2B_ cos 2uJot - 2^2 sin^ ujot 



-1, 
1 



2e I C_ cos ujQt I 
which indicate that when t ^ (^"-i)?^ ^^ = 1,2,...) and 



(39) 
(40) 



2uio 
2 



1 — 2B± cos 2LjQt — 2C± sin u^t < 2e I C± costuot 



(41) 



we have Sj{t) < and Sj{t) < which means that the supercurrent exhibits squeezing. 
However, when t = ^^'^^^-'^ (n = 1,2,...), we have Sj{t) = Sj{t) = +oo. From the 
uncertainty relation ( pTf ) and Eq. (^) we know that 



((A/)^)((Ag)2)>0 



(42) 
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which means that no squeezing of the supercurrent occurs under this circumstance. 

We plot the time evolution of squeezing of the supercurrent in the MJJ when 
cuq — 2NuJi for even and odd CS cases in Fig.l and Fig. 2. As can be seen from Fig.l, 
in the case of even CS, the squeezing is weakened with the increase of the value of 
\Z\ for a given A^. Nevertheless, the case of odd CS is contrary to the even CS case. 
From Fig. 2 we see that the squeezing is strengthened with the increase of the value 
of I Z I for a given N , although there does not exist squeezing of the supercurrent 
when Z is very small. It is noted that the singular points of the squeezing parameters 
appear periodically with the period ^ which depends only on the external classical 
voltage bias, the number of the singular points increases with the increasing of N . 
Therefore, the angular frequency of the nonclassical EMF does not affect periodicity 
of the supercurrent squeezing. Comparing Fig.l with Fig. 2, we find that the squeezing 
curve of the odd CS approaches to that of the even CS with the increasing of |Z | for 
a given N. 

Since even and odd CS are nonclassical states while the Glauber CS are classical 
states, their effects on the squeezing of the supercurrent in the MJJ are quite different. 
To see this, we write out the squeezing degree of the supercurrent for Glauber CS, 



Si{t) = 



(e2^^ — e 2^^){1 + e ' 


cos[Ai \Z 


sin(a'it — < 




2e 


cos [2^ 


\Z\ 


sin(a;i - 


- Lp) + UqI] 





- 1 (43) 



which indicates that when 2^ \Z \ sm.{ijj]t — (p) + ujot ^ {2n — l)7r/2 and (e^^^ — e~^^^) 
{1 + e~^^ cos[4^ \Z I sm{uJit — ip) + 2uJot]} < 2e | cos[2^ \Z \ sin(u;i — 99) + uiQt] \ the 
supercurrent exhibits squeezing (i.e., Si{t) < 0). When 2$, \Z \ sin(a;i — if) + cuot — 
(2n-i)7r ^^ _ ]^ 2, . . .), the squeezing disappears due to Sj{t) — +00 
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FIG.l: Time evolution of squeezing of the supercurrent in the MJJ 
when the external nonclassical EMF is in even CS for </? = ^ and 

(a) \Z\= 0.5, iV = 5; (b) \Z\= l,N = f>- (c) \Z\= 0.5 iV = 1; (d) \Z\= 0.1 N = 1. 
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FIG. 2: Time evolution of squeezing of the supercurrent in the MJJ 
when the external nonclassical EMF is in odd CS for (/9 = | and 

(a) \Z\= 0.5, iV = 5; (b) \Z\= l,N = f>- (c) \Z\= 0.5 iV = 1; (d) \Z\= 0.1 N = 1. 
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FIG. 3: Time evolution of squeezing of the supercurrent in the MJJ 
when the external quantum EMF is in the Glauber CS for = ^ and 

(a) \Z\= 10, = 5; (b) \Z\= 10, iV = 1; (c) |Z|= 1 = 1. 

In Fig. 3, we plot the time evolution of the squeezing of the supercurrent in one 
period in the presence of the Glauber coherent-state EMF when the step condition is 
satisfied, i.e., ujq = Nui. Obviously, the squeezing of this case is quite different from 
that of even and odd CS. From (a) and (b) in Fig. 3, we can see that for a given \Z \ , 
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oscillations of Si{t) are strenghtened during the preceding half period. From (b) and 
(c) we can see that for a given value of A^, the oscillatary behaviors are weakened and 
the number of the singular points in every period decreases with the decreasing of |Z |. 

4 Quantum coherences of supercurrent in the MJJ 
in the presence of quantum EMF 

In this section we consider another important aspect of statistical properties of the 
supercurrent in the MJJ, that is its quantum coherences which describe the degree 
of quantum phase correlation in the supercurrent. The degree of the second-order 
quantum coherence can be obtained by evaluating the second-order correlation function 
of the supercurrent at two different spacetime points. If purely temporal correlations 
are of interest, the relevant correlation function in the normalized form can be defined 
by 

)_W<±lM (44) 

{/(*)){/(* + T)) 

where the bracket "( )" denotes taking expectation values over states of the external 
quantum EMF, and ": :" stands for normal ordering of quantum mechanical operators. 
Following considerations in quantum optics. We are interested in behaviors of the 
second-order correlation function (Kl) at r = 0. Then, from Eq. (PI) we get 



,(^)(t)^,(^)(M) = ii||l^ (45) 

Substituting the expression of the supercurrent operator (^ into the above equation 
and making use of the formula 

: D{2^e''''') := e^^' D{2^e''^'') (46) 

we find the general expression of the second-order correlation function in the following 
form 
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For even and odd CS, we find that 

1 — 2i?_|_ cos 2ti;oi (2) 1 — 2i?„ cos 2ci;o^ 
2C+ sm ujot 2C_sin ujqi 

where B± and C± are given by Eq. (PD|). 

From Eq. (^), it can be seen that when t = ^{n = 0, 1, 2, . . .) we have gi^\t) = 

g^\t) = +00 which indicates that quantum coherences vanish since the mean value of 

the supercurrent is zero at these points. 



FIG. 4: Time evolution of quantum coherence of the supercurrent 
when the external nonclassical EMF is in even CS for 9? = § 

and (a) \Z\= 0.5, N = 1; (b) \Z\= 0.5, N = 2. 

In Fig. 4 and Fig. 5 we plot the time evolution of quantum coherences of the super- 
current in the MJJ when uq = 2Nuji for even CS and odd CS, respectively. As can be 
seen from these figures that when the external nonclassical EMF is in even and odd 
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CS, the supercurrent in the MJJ exhibits similar quantum coherences during the time 
evolution. 



FIG. 5: Time evolution of quantum coherence of the supercurrent 
when the external nonclassical EMF is in odd CS for = ^ 

and (a) \Z\= 0.5, = 1; (b) \Z\= 0.5, = 2. 

It is worthwhile mentioning that when the external quantum EMF is in Glauber 
CS, the time evolution of quantum coherences of the supercurrent in the MJJ is quite 
different from that of the even and odd CS cases. For the case of the Glauber CS, 
it is quite easy to find the second order correlation function to be g^'^\t) = = 1 
which means that the quantum coherence of the supercurrent in the presence of the 
coherent-state EMF does not change during the time evolution. Therefore, nonclassical 
states and classical states of the external quantum EMF lead to quite different time 
evolution of quantum coherences of the supercurrent in the MJJ, in this sense, the MJJ 
can distinguish the external nonclassical hght fields and classical light fields. 
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5 Concluding remarks 

We have studied some properties of the supercurrent in the MJJ in the presence of 
nonclassical hght fields, the CVSS, squeezing effect and quantum coherences. In par- 
ticular, we have investigated in detail these properties when the external quantum EMF 
is in even and odd CS. We have found that the CVSS of the MJJ depends strongly on 
the off-diagonal elements of the density matrix which describes the external quantum 
EMF. We have shown that the supercurrent in the MJJ exhibits both squeezing and 
quantum coherences. However, when the CVSS of the MJJ is satisfied, squeezing and 
quantum coherences vanish periodically due to the periodical appearance of singular 
points of the squeezing parameter and the second-order correlation function of the su- 
percurrent. we have also shown that the MJJ can feel the difference not only between 
classical light fields and nonclassical light fields but also betwwen different nonclassical 
light fields. This indicates that the MJJ are very sensitive devices that can respond 
to the quantum nature of external nonclassical hght fields. Therefore, the MJJ is a 
sensitive detector to nonclassical hght fields. We believe that the results obtained in 
this paper could be useful in the study of the controlling of Josephson tunnelling and 
the detection of nonclassical light fields in quantum optics. 

The research was supported in part by the National Nature Science Foundation of 
China. 
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